Bundling up carbon nanotubes through Wigner defects 
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We show, using ab initio total energy density functional theory, that the so-called Wigner defects, 
an interstitial carbon atom right besides a vacancy, which are present in irradiated graphite can also 
exist in bundles of carbon nanotubes. Due to the geometrical structure of a nanotube, however, this 
defect has a rather low formation energy, lower than the vacancy itself, suggesting that it may be 
one of the most important defects that are created after electron or ion irradiation. Moreover, they 
form a strong link between the nanotubes in bundles, increasing their shear modulus by a sizeable 
amount, clearly indicating its importance for the mechanical properties of nanotube bundles. 



Single-walled carbon nanotubes (SWNT) have at- 
tracted an enormous amount of attention since its 
discovery^ about ten years ago. No doubt this is re- 
lated to all their extraordinaire materials properties and 
potential for applications 2] . From a mechanical point of 
view, SWNT have exceedingly large elastic modulus Q 
around 1 TPa, which would make them perfect for low- 
density high-modulus fibers. When in bundles, however, 
a much worse performance is usually obtained, mostly 
related to the weak interaction between tubes, which 
makes it relatively easy for them to slide against each 
other£|. A breakthrough has been recently obtained^ 
using electron-beam irradiation to generate crosslinks be- 
tween the tubes. This increased the shear modulus by a 
factor of approximately 30 thus indicating a new way to 
fabricate stronger fibers. 

Irradiation by electrons or ions is being used by many 
groups [3 as an important tool to modify the structure 
and properties of nanotubes. As examples, we can men- 
tion the welding of tubes 6], and the dramatic increase 
of shear modulus for a bundle of nanotubes In or- 
der to have a better control of the nanotube engineering, 
it is fundamental to have a detailed microscopic under- 
standing of the defects that are created upon irradiation. 
Properties like formation energy, geometry, stability and 
electronic structure are critical in determining what type 
of defects will be produced and how they will affect the 
mechanical behavior of the carbon fiber. Some theoreti- 
cal works have investigated the effects of electron or ion 
irradiation on SWNT and SWNT bundles @, 0. These 
studies consisted of Molecular Dynamics simulations with 
the interaction between the atoms being described by em- 
pirical potentials. Even though these works are relevant 
and indicative of the processes that may occur, they may 
miss important aspects since empirical potentials cannot 
describe appropriately features such as re-hybridization, 
rebonding or charge redistribution, all of them quite im- 
portant when chemical bonds are being broken and made, 
specially in stressed configurations that may result from 
the irradiation procedure. Moreover, recombination bar- 
riers and formation energies are usually poorly described 
by these potentials. Yet another difficulty with these 
potentials is the important fact that carbon atoms have 
many possible hybridizations. Therefore, even though 



cross links between tubes in bundles have already been 
predicted to appear after irradiation |9j, their detailed 
structure and properties that will control how they affect 
the mechanical behavior of bundles have not yet been de- 
termined. It is thus clear that a fully quantum mechani- 
cal description of the electronic structure of the system is 
necessary, and that is precisely what we describe below. 

All our results are based on ab initio total energy den- 
sity functional theory^fj calculations. We have used 
ultrasoft pseudopotentials|llj|. a plane wave expansion 
up to 230 eV, and the local density approximation for 
the exchan ge-c orrelation potential as implemented in the 
VASP code 12]. A defect-in-supcrcell approach was em- 
ployed, similar to what has been used in the investigation 
of defects in graphite[l3L Il4). We have in our supercell 
two (5,5) nanotubes, with 100 atoms each, with lattice 
primitive vectors appropriate for a hexagonal arrange- 
ment of tubes in the bundle, as shown in Fig. 1(a). 
The Brillouin zone sampling had a total of six points, 
with three of them along the tube axis direction. In all 
calculations the positions of the unconstrained atoms in 
the supercell were relaxed until the forces in them were 
smaller than 0.02 eV/A. We obtain 1.42 A for the C- 
C bonds perpendicular to the tube axis, and 1.41 A for 
the others, which gives for the cell length the value of 
12.21 A. For the inter-tube distance in the bundle we 
obtain approximately a center-to-center distance of 9.87 
A. From a variety of tests we estimate that the errors in 
energy differences due to all our approximations are of 
the order of 0.05 eV. 

The two tubes in the unit cell are arranged in such 
a way as if they were obtained by rolling up two neigh- 
boring graphene sheets in a graphite structure, in other 
words, at the "line-of-contact" between them the atomic 
arrangement resembles that of graphite (see Fig. 1(b)). 
Since the inter-tube interaction is weak, they have a 
lot of rotational freedom within the bundle at room 
temperature[15j, which does not make this choice restric- 
tive. However, the curvature will play an important role 
in the stability of the defect, as discussed below. The 
structures of a variety of defects in graphite have been re- 
cently studied |13lll4| . In particular, an intimate Frenkel 
pair defect has been suggested to be responsible for the 
undesirable energy release in graphite moderators in nu- 
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FIG. 1: (a) Final geometry for a SWNT bundle, showing the 
hexagonal order and the unit cell (tubes in black), (b) Detail 
of the "line-of-contact" between the tubes in the defect free 
configuration (region marked by the rectangle in (a)), as seen 
from the inside of the top tube. The orange (grey) atoms are 
in the top (bottom) tube. The atom that will form a link 
between the two tubes is marked in red, whereas its three 
nearest neighbours are marked in blue. Upon displacement, 
the red atom will form bonds to the two atoms in the top 
tube marked in green, (c) Same atoms after the formation 
of the defect, (d) Final geometry for the configuration with 
the Wigner defect. The interstitial atom (red) can be clearly 
seen connecting two tubes, (e) Detailed view of the Wigner 
defect. The pentagon formed after the interstitial is ejected 
can be clearly seen (atoms 7 through II). 



"line-of-contact" and curvature were both favorable. 




FIG. 2: Charge density (e~/A ) in: (a) plane perpendicu- 
lar to the tubes axis and farther away from the Wigner de- 
fect. Five C-C bonds in each tube can be clearly seen, (b) 
Plane passing through atoms 1-2-3-4 in Fig. 1(e). The C-C 
bond that existed between atoms 1 and 2 in the upper tube 
is broken, and two new bonds are formed with the interstitial 
atom, (c) Plane that approximately passes through the five 
atoms marked in Fig. 1(e) from 7 to 11. (d) Similar to (c) in 
graphite. 



clear reactors, the Wigner energy. The procedure we em- 
ployed in the attempt to produce a Frenkel pair similar 
to what was found in graphite |l3tll4j |. was the following: 
we moved an atom in one of the tubes from the "line- 
of-contact" and placed it approximately in between the 
two tubes, in the middle of a close by C-C bond in the 
other tube. The C atoms in the farthest away planes (in 
both tubes, 20 atoms total) from the defect were held 
fixed throughout the atomic relaxations. This was done 
in order to simulate an isolated defect, where the atoms 
far from the perturbation would not move. During the 
first few relaxation steps the displaced atom was also 
kept fixed, otherwise it would tend to return to its orig- 
inal position. All the attempts to create a Frenkel pair 
defect by displacing one atom from the top tube of Fig. 
1(b) (one of the orange atoms) were unsuccessful. Due 
to the curvature of the tube, the atom always returned 
to its original position. Eventually, tube rotations may 
be able to stabilize such configurations, but this may de- 
pend on the temperature and the time scales for rotation 
as compared to the defect relaxation time. On the other 
hand, when we displaced an atom from the bottom tube 
of Fig. 1(b), a stable Frenkel pair was obtained since the 



As shown in Fig. 1, the displaced atom (marked in red 
and labeled 3 in Fig. 1(e)) forms two bonds with C atoms 
of the top tube (marked in green, and labeled 1 and 2 in 
Fig. 1(e)). The bond length between atoms 1 and 2 in- 
creases from their equilibrium value of 1.42 A to 2.16 A, 
which clearly indicates that it is broken. This can also 
be seen in Fig. 2 via the charge density plots. On the 
other hand, the distances between carbons 1-3 and 2-3 
(Table I) clearly show the formation of two C-C bonds. 
Curvature effects and the constraints being imposed lead 
to different bond lengths. The displaced atom remains 
bonded to its original tube through one of its original 
neighbors, the atom labeled 4 in Fig. 1(e). The distance 
3-4 decreases from 1.41 A to 1.36 A upon the formation 
of the defect. This is suggestive of a sp 2 hybridization 
for these two carbon atoms with the formation of a dou- 
ble bond between them, since the C-C distance in the 
H2C=CH2 molecule is 1.34 A. These results indicate that 
a strong bond is being formed between the two tubes, as 
can be seen from the charge density plots of Fig. 2. The 
atom 4 remains linked to its two other neighbors via the 
bonds 4-5 and 4-6, with bond lengths diplayed in Table 
I. As a reference for comparison, we have also studied the 
Wigner defect in graphite [lall4| . We have used a similar 
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TABLE I: Interatomic distances (A) for the atoms labelled 1 
through 11 in Fig. 1(e). Results are presented for the SWNT 
bundle when some of the atoms were fixed (line I) , when they 
were all allowed to relax (line II), and for graphite (line III). 
1-2 1-3 2-3 3-4 4-5 4-6 7-8 8-9 9-10 10-11 7-11 

I 2.16 1.50 1.41 1.36 1.42 1.50 1.57 1.41 1.41 1.40 1.42 

II 2.21 1.45 1.43 1.33 1.44 1.45 1.54 1.42 1.40 1.40 1.43 

III 1.57 1.45 1.44 1.32 1.44 1.44 2.00 1.39 1.40 1.40 1.39 



approximation as described above, however the supercell 
now had 100 C atoms, 50 in each of the graphite planes, 
and 100 k-points were used for the Brillouin zone sam- 
pling. The equilibrium C-C bond length was 1.41 A and 
the interplanar distance was 3.30 A. Upon the formation 
of the Frenkel pair, we obtained the equivalent distances 
shown in Table I, which indicate that the links caused by 
the defects in the nanotube bundles will most likely be 
stronger than in graphite. 

An important aspect of this defect in nanotube bun- 
dles when compared to graphite is the atomic relaxation 
around the defect. When the C atom is displaced, it re- 
mains bound to one of its three former nearest neighbors 
(the atoms marked in blue in Fig. 1), as discussed above. 
The other two C atoms are now in an energetically un- 
favorable configuration, and will tend to approach each 
other in an attempt to make a chemical bond. This will 
lead to the formation of a pentagon, which can be clearly 
seen in Fig. 1(e), where the five C atoms are labeled 
from 7 through 11. In a nanotube, however, due to its 
curvature allied to a certain degree of flexibility towards 
compression, these two C atoms can get relatively close 
to each other, with a final interatomic distance of 1.57 
A. This should be contrasted to the similar distance in 
graphite, 2.00 A. Therefore, in the case of nanotubes 
these two carbon atoms can rebond much more effec- 
tively, which should have an important contribution to 
the stability of the defect. In Fig. 2(c) and Fig. 2(d) we 
present charge density plots in planes that pass through 
the atoms that form this pentagon, for both the nan- 
otube and graphite. It can be clearly seen that there 
is indeed a much stronger rebonding in the nanotubes 
when compared to graphite. It should be noted that due 
to the curvature of the nanotube, it is very difficult to 
have a plane passing exactly through all five C atoms 
(from 7 through 11), that is the reason why atom 9 (see 
Fig. 1(e)) appears less clear in Fig. 2(c). All the relevant 
interatomic distances are summarized in Table I. 

We calculate the defect formation energy as the differ- 
ence between the total energy for the system with and 
without the defect, which gives a value of 7.27 eV. For 
the Wigner defect in graphite we obtain a formation en- 
ergy of 10.67 eV, similarly to other authors [HI ll^ . This 
means that it is much easier to produce these defects 
in nanotube bundles than in graphite, indicating that it 
may be very common in irradiated samples. If we con- 
sider a vacancy in an isolated nanotube, we obtain for 



a fully relaxed structure a formation energy of 5.82 eV. 
However, if we now also allow all the atoms in the bundle 
to relax, i.e., the 20 C atoms that were held fixed can now 
relax, we obtain a formation energy of 5.53 eV. This is 
an extremely low value, showing that this defect may be 
one of the most important defects in nanotube bundles. 
The full relaxation allows the tubes to orient themselves 
in such a way that their interlink bonds are stronger, as 
indicated by the bond lengths shown in Table 1. There- 
fore, if the tubes have enough mobility to accommodate 
the defect, or if the defect concentration is high enough, 
we expect a formation energy close to 5.5 eV. However, 
if the defect is isolated and the tubes cannot freely move 
the formation energy may be closer to 7.3 eV. 




Displacement (A) 

FIG. 3: Change in energy (per unit length) caused by the 
relative displacement of two isolated tubes (not in a bundle). 

The strong link that is being formed between the tubes 
may have a dramatic change in the mechanical proper- 
ties of nanotube bundles. To estimate the changes in 
the shear modulus, we perform calculations for two iso- 
lated tubes, instead of a full bundle. This should give 
the correct trend, since the tube-tube interaction is much 
stronger in the presence of the Wigner defect. We per- 
form relative displacements as illustrated in Fig. 3. For 
the displacements perpendicular to an imaginary line 
that joins the tubes axes we obtain the dashed curves, 
whereas for displacements making a 30° angle with re- 
spect to the previous displacement direction we have the 
solid curves. This latter direction would be more rep- 
resentative of a shear displacement of tubes in a bun- 
dle. For the dashed curves we obtain shear moduli of 3 
GPa and 59 GPa for the pure and defected tubes, respec- 
tively, whereas for the solid curves we obtain 7 GPa and 
97 GPa, respectively. This indicates an increase in the 
shear modulus by a factor of 20 in the former case and by 
a factor of 14 in the latter case. Therefore, we show that 
this defect may be the main reason behind the observed 
increase in shear modulus of nanotube bundles upon 
irradiation 0. In the work by Kis and collaborators 
it has been proposed that carbon interstitials or radia- 
tion induced chemical reactions involving, for example, 
carboxyl groups could form links between the nanotubes 
in a bundle. These changes would be responsible for the 
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increase in shear modulus, however, no estimates for the 
formation energies of these structures are provided. We 
here propose a similar idea but with an alternative local 
atomic structure, which due to its low formation energy 
is very likely to be created during the irradiation. 
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FIG. 4: Energy as a function of displacement of the intersti- 
tial atom from its ideal configuration (zero displacement). 

Finally, an important question is related to the stabil- 
ity of such a defect. In graphite, it has been proposed [l3j 
that this defect is responsible for the Wigner energy re- 
lease peak at 200°C, associated with an experimental 
barrier^(| of approximately 1.4 ± 0.4 eV, which agrees 
nicely with the calculated barrier [T3HT4] ] of 1.3 eV for the 
defect recombination. To estimate the barrier, we lin- 
early interpolate the coordinates of the displaced atom 
between its position in the pure bundle (Fig. 1(a)) and 
in the final defect configuration (Fig. 1(d)). We then 
held this atom fixed and let all the other atoms relax, 
except the 20 constrained atoms as discussed above. Ten 
intermediate configurations were calculated, and the re- 



sults are plotted in Fig. 4. We obtain a recombination 
barrier of approximately 2.4 eV, indicating that the de- 
fect is significantly more stable in nanotube bundles than 
in graphite. If we assume that the recombination of this 
defect would also be associated with an energy release, 
like in graphite, and that to observe such an energy re- 
lease it would require a recombination rate of the same 
order than that in graphite, we estimate that in nanotube 
bundles such a release should be observed around 600°C. 

In summary, we have identified a very stable Frenkel 
pair in nanotube bundles, able to provide strong links 
between nanotubes, thus altering in a dramatic way the 
mechanical properties of these bundles. Since these de- 
fects have rather low formation energies, they should be 
prevalent in irradiated samples. They should play an 
important role in nanotube engineering using irradiation 
sources, a field that has been gaining a lot of attention 
lately. Besides forming the strong links, as already men- 
tioned, they could be spots where welding of two tubes 
will start. Moreover, since they provide a somewhat open 
local structure in one of the tubes (a small "hole" in the 
tube wall) , they may allow the entrance of atoms or small 
molecules inside the tube, which may be important for 
storage purposes 0. Even though we focused on a (5,5) 
nanotube, the conclusions are not dependent on its chi- 
rality. Finally, it has not escaped our attention that a 
similar structure will probably exist in multi- walled nan- 
otubes linking, for example, two neighboring walls |Tsj. 
Some of the ideas discussed here were already foreseen in 
the study of defects in graphite^], however, the struc- 
ture of nanotubes is such that the defect has a combined 
high stability and low formation energy, which could not 
be expected from graphite studies. 
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